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Abstract 

In this paper, using a novel approach involving the truncated Laurent expansion 
in the Painleve analysis of the (2+1) dimensional K-dV equation, we have trilin- 
earized the evolution equation and obtained rather general classes of solutions in 
terms of arbitrary functions. The highlight of this method is that it allows us to con- 
struct generalized periodic structures corresponding to different manifolds in terms 
of Jacobian elliptic functions, and the exponentially decaying dromions turn out 
to be special cases of these solutions. We have also constructed multi-elliptic func- 
tion solutions and multi-dromions and analysed their interactions. The analysis is 
also extended to the case of generalized Nizhnik-Novikov-Veselov (NNV) equation, 
which is also trilinearized and general class of solutions obtained. 



1 Introduction 



The recent interest in the investigation of integrable models in (2+1) di- 
mensions can be attributed to the identification of dromions in the Davey- 
Stewartson I equation [TJ[2] . These dromions which originate at the crosspoint 
of the intersection of two nonparallel ghost solitons decay exponentially in 
all directions and are driven by lower dimensional boundaries or velocity po- 
tentials (arbitrary functions) in the system. It must be mentioned that the 
presence of these so called boundaries enriches the structure of (2+1) dimen- 
sional integrable models leading to the formation of localized solutions like 
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dromions, lumps, breathers, positons, etc. [2|3] . Thus, a judicious harnessing 
of these lower dimensional arbitrary functions [4,5 6 7 8,9] of space and time 
may give rise to new spatio-temporal patterns in higher dimensions besides 
throwing more light on their integrability. 



In this direction, we have recently identified a simple procedure to solve the 
(2+1) dimensional systems, namely the Painleve truncation approach (PTA) 
[7]. This approach converts the given evolution equation into a multilinear (in 
general) equation in terms of the noncharacteristic manifold. This multilin- 
ear equation can be solved in terms of lower dimensional arbitrary functions 
of space and time. Using this approach, one can generate various generalized 
periodic and localized solutions. This was earlier demonstrated for the case 
of long wave-short wave resonance interaction equation in (2+1) dimensions 
by converting it into a trilinear equation [7]. In this paper, we wish to show 
that through the PTA the (2+1) dimensional K-dV equation equation can 
be trilinearized. A four parameter special class of solutions involving three 
arbitrary functions g(y), h(y) and f(x,t) in the indicated variables naturally 
arises for this trilinear equation. This in turn leads to the universal form of 
solution obtained earlier by Tang, Lou and Zhang [5] using a variable separa- 
tion approach as a special case. However, there exists a possibility of obtaining 
more general solution to the trilinear equation, which remains to be explored. 
A large class of Jacobi elliptic function periodic solutions, multidromion and 
bound state solutions are also obtained. It is also shown that the procedure 
can be extended to the case of Nizhnik-Novikov-Veselov (NNV) equation as 
well. 



The plan of the paper is as follows. In section 2, we show how the (2+1) di- 
mensional K-dV equation can be trilinearized using the PTA and the solution 
obtained in terms of three arbitrary functions and four arbitrary parameters. 
Sections 3 and 4 contain two broad classes of solutions of the (2+1) dimen- 
sional K-dV equation, both periodic and exponentially localized ones, through 
judicious choices of the arbitrary functions. In section 5, we point out how the 
generalized NNV equation can be trilinearized through the PTA and how var- 
ious periodic and localized solutions can be obtained. Finally in section 6, 
we summarize our results. The appendix contains the one dromion solution of 
the (2+1) dimensional K-dV equation obtained through Hirota bilinearization 
approach for comparison. 
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2 The (2+1) dimensional K-dV equation and construction of solu- 
tions 



The (2+1) dimensional K-dV equation introduced by Boiti et al. [10] has the 
form 

Qt + Qxxx = ^(qd~ 1 q x ) x . (1) 
This nonlocal equation (jTJ) reduces to the K-dV equation if x = y. Introducing 
a potential function v(x,y,t) defined by the equation 

v y = q x , (2) 
Eq. ([T]) can be converted into a set of coupled equations of the form 



qt + q X xx = 3(qv) x , (3a) 

v y = q x . (3b) 

The nature of f!3bl) permits the presence of an arbitrary function v'(x,t) as 

ry 

v(x,y,t) = / q x dy+v(x,t). (4) 

J — oo 

Expanding the physical field q and the potential v in the form of a Laurent 
series in the neighbourhood of the noncharacteristic manifold <p(x,y,t) = 0, 
(<f) x = <p y 0) and utilising the results of the Painleve analysis [TT], we obtain 
the following Backlund transformation by truncating the Laurent series at the 
constant level term 



q = q (j) 2 + qi4> 1 + q 2 , (5a) 
v = v (fi~ 2 + vi^ 1 + t> 2 , (5b) 

where (q,v) and (52,1^2) are two different sets of solutions of the (2+1) di- 
mensional K-dV equation. In ref. [1] certain basic localized dromion solutions 
have been obtained using (jSJ) and Hirota bilinearization. Here, we make use 
of the full power of Eq. (j5J) by treating it as a variable transformation, and 
obtain rather general classes of solutions. This new approach leads to a wide 
class of solutions not only for Eq. (pQ) but for other (2+1) dimensional systems 
as well [7]. 

Considering now a vacuum solution of the form 

q 2 = 0,v 2 = v 2 (x,t) (6) 

and substituting (jSJ) with (jSJ) in Eq. we obtain the following set of equa- 
tions by equating the coefficients of (0 -5 , 0~ 3 ) to zero, 
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- 24q <f>l + 12q v <l) x = 0, (7a) 
-Vo4> y + q (fix = 0. (7b) 

Solving the above system of equations, we obtain 



q = 2cj) x (j)y, (8a) 

= 24> 2 X . (8b) 

Again, collecting the coefficients of (0~ 4 , 0~ 2 ), we have 

18qox(j>l + ^8q (j) x (f) xx - 6gi0^ - 3(q 0x v - 2q vi4> x 

-qi<PxV ) - 3(v 0x q - 2v q 1 (p x - v^qo) = 0, (9a) 

v 0y - v x <\) y - (q 0x - q X (j) x ) = 0. (9b) 

Solving the above system of equations, we obtain 

qi = -2(f) xy , (10a) 

«i = -2(j) xx . (10b) 



Next, collecting the coefficients of (0 3 , x ), we have 



- 2q (f) t + {-Qq 0xx (p x - Gqo x (fixx ~ 2q (j) xxx + 6q x 
+6qi4> x (t>xx) ~ 3(q 0x V! - 2q (f> x v 2 + q lx v - qiVx(f> x ) 

-3(qoVi x + qiv 0x - qiV X (j) x ) = 0, (11a) 
vi y -qi x = 0. (lib) 



Here (lllbj) is an identity as may be verified using (jTUj). Solving filial) by using 



Yt yxxx , Yxyxxv YxxYxy tic\\ 

v 2 = — tt, 1 , — ; • ( 12 J 



and (ITU1) . we get the form of v 2 as 

+ 

30a; > n 

Now, we collect the coefficients of next order (0~ 2 , 0°) to obtain 

(go* - gi0i) + (qoxxx - 3qix4>xx - qi4>xxx - 3qixx4> x ) 

-3(qoxV 2 + qixVi - qiv 2 (j)x) ~ 3(g t> 2 z + qiVi x ) = 0, (13a) 

v 2y = 0. (13b) 



4 



On the other hand, substituting the expression for the quantities go, vq, Qi, v± 
and v 2 in Eq. (j 13 aft , we obtain the trilinear equation 

$xx(jPy&xyy 4 > yy4 > xy) ( f ) x( < fiyy4 ) xxy 4*y4*xxyy) 0- 

(14) 

Substituting ffT2l) in (I13bl) . we obtain 



(<Pty + <j>xxxy)<l>x<t>l ~ (4>t + (f>xxx) (f>xy4>l + 3((j) x (j) XX yy - <p xx (j) X yy)<p x (i) y 

-3((j) x (f) xx y - (j) xx (j) x y)(f) y (j) x y - 3((j) x (f) xx y - (f) xx (f) x y)(j) x (f)y y = 0. (15) 

Using Eq. (|T4"|) into (|T5|) . the later reduces to the form 



(4>ty + 4>xxxy)4>x4>y ~ (4>t + 4>xxx) 4>xy4>y ~ 3{4>x4>xxy ~ <Pxx(f>xy)4>xy = 0. (16) 

which is again in a trilinear form. Thus equations (fl4"|) and (|T6|) may be con- 
sidered as the equivalent trilinear forms of Eqs. flS}. 

One can immediately observe that the set of trilinear equations (|14p and ffTBT) 
admit a set of two arbitrary functions and their products as solutions 



(j> = f(x,t), (17a) 

= g(y) (17b) 

and 

<j) = f{x,t)g{ y ), (17c) 

where / and g are arbitrary functions in the indicated variables. In fact a 
more general solution involving three arbitrary functions and four arbitrary 
constant parameters can be easily identified: 

<f)(x, y, t) = a + c 2 f(x, t) + c 3 h(y) + c 4 /(x, t)g(y), (18) 

where h(y) and g(y) are in general different arbitrary functions of y. Here ci, 
C2, C3, and C4 are arbitrary parameters. The problem of identifying even more 
general solution to (j!4j) and (TT6"j) remains to be investigated, while we will 
concentrate on the special solution f|T8|) in this paper. 



Collecting the coefficients (0 we obtain 

Qit + Qixxx ~ 3qixV 2 ~ 3qiv 2x = 0. (19) 

It can be checked that equation f|T9l) is compatible with the earlier results. 

Now substituting the above form (|T8l) for the manifold <p(x, y, t) into the trun- 
cated Painleve series §5§ for the functions q(x,y,t) and v(x,y,t), along with 



5 



the expressions for the coefficient functions qo, vq, qi, V\ and v 2 given above, 
we finally obtain the solution to Eqs. ([3]) as 



/ ,n = Vx[{c2 + c^g)c 3 h y - (ci + c 3 fe)c 4 ff y ] 
^' V ' } ~ h + c 2 f(x, t) + c 3 %) + c 4 /(x, ^(y)]a ' 1 J 

v(x,y,t)- 



[ci + ca/(x, t) + c 3 %) + c 4 /(x, t)g{y)\ 
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2(C 2 + 04(7)/^ + ft + fxxx ^ 



c\ + c 2 f(x, t) + c 3 h(y) + c 4 f(x, t)g(y) 3f x 



It may be noted that the expressions given in (|20p and (I2T|) coincide exactly 
with the forms obtained by Tang, Lou and Zhang [5] using the method of 
variable separation for the special case g(y) = h(y) in Eq. ([181) . Another special 
case C\ = c 2 = 0, c 3 = c 4 = 1 was identified by Peng recently [9]. Also, in our 
case it is clear that finding any solution to (1141) and (fl6l) which is more general 
than (flSI) will lead to more general solution than (l20l) -( l2TT) . This problem is 
being investigated further. In the following, we will investigate the nature of 
solutions (T20l ) - fT2TT) corresponding to periodic and localized solutions. For this 
purpose, we shall consider the special cases (i) c 4 = and (ii) c 2 = c 3 = in 
(fT8l) or (j20|) - (j2ip . corresponding to the sum and product of arbitrary functions, 
respectively, and investigate the nature of periodic and localized solutions. It 
may be noted that there is no difficulty in proceeding with the general form 
(!20l) - (pT|) also except for the fact that the expression will be more lengthy 
which we desist from presenting here. 



3 Periodic and localized solutions corresponding to sum of arbi- 
trary functions 



In this section, we will assume c\ = c, c 2 = c 3 = 1, c 4 = in <HM or (j20|) - (j2T|) . 
Consequently, we have 



^ f x ^fxx 1 ft ~H fxxx (22b) 



U + h + cf (f + h + c) 3f x 
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3.1 Harnessing of arbitrary functions and novel solutions of (2+1) dimen- 
sional K-dV equation 



Let us now choose the arbitrary functions f(x,t) and h(y) to be Jacobian el- 
liptic functions, namely sn or cn functions. The motivation behind this choice 
of arbitrary function stems from the fact that the limiting forms of these func- 
tions happen to be localized functions. Hence, a choice of cn and sn functions 
can yield periodic solutions which are more generalized than exponentially 
localized solutions (dromions). We choose, for example, 



f = a sn(aa; + C\t + d\\ mi), h = (3 sa(by + d 2 ; m 2 ) (23) 



so that 



2a P abcn(ui, m 1 )dn(u 1 ; 7%) cn(-u 2 ; m 2 )dn(u 2 ; m 2 ) ,_, N 

q{x, y, t) = — vnr~t > ( 24 ) 

(c + asn{ui; mi) + (3sn{u 2 ; m>2)) 



where u% = ax + c\t + d\ and u 2 = by + d 2 . In Eqs. ( 1231) and f l24l) . the quanti- 
ties mi and m 2 are moduli of the respective Jacobian elliptic functions while 
a, (3, a, b, c, ci, d\ and d 2 are arbitrary constants. We choose the constant 
parameters c, a and (3 such that |c| > \a + j3\ for nonsingular solution. The 
profile of the above solution for the parametric choice a = (3=1, a=0.5, 6=0.4, 
c=-4, Ci=-1, <ii=(i 2 =0, m!= 0.2, m 2 = 0.3 at time t=0 is shown in Fig.l(a). 
Note that the periodic wave travels along the x-direction only. 



3.1.1 (1,1) dromion solution 

As a limiting case of the periodic solution given by Eq. (1241) . when mi, m 2 
— > 1, the above solution degenerates into an exponentially localized solution 
(dromion). Noting that cn(u; 1) = dn(u; 1) = sechu and sn(u; 1) = tanhu, the 
limiting form corresponding to (1,1) dromion takes the expression 

2a b a (3 sech 2 (ax + c\t + di)sech 2 (by + d 2 ) (| | > | + (3\) (25) 



(c + atanh(ax + c\t + di) + (3tanh(by + d 2 )) 2 ' 
and the variable v then takes the form (using expression f!22bl) ) 



2a 2 a 2 sech (ax + c\t + d\ 



(c + atanh(ax + C\t + d\) + /3tanh(62/ + d 2 )) 2 
4a a 2 sech 2 (ax + c\t + di)tanh(ax + cit + d\) 
(c + a tanh(ax + C\t + d\) + (3t&nh.(by + d 2 )) 



+ 



— [ci + 4a 3 tanh 2 (ax + c x t + di) - 2a 3 sech 2 (ax + c x t + d x )]. (26) 

oCL 
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Schematic form of the (1,1) dromion for the parametric choice a = (3=1, 
a=0.5, 6=0.4, c=-4, Ci=-1, di=d 2 =0 at time t=0 is shown in Fig. 1(b). Again 
note that the dromion travels along the x-direction. One can obtain the (1,1) 
dromion given by Radha and Lakshmanan in ref. jl] by fixing the parameters 
a, (3, a, b, c\, d\, d 2 of equation ([25]) suitably (see also the Appendix). 



3.1.2 More general periodic solution and (2,1) dromion 

Next we obtain a more general periodic solution by choosing further general 
forms for the arbitrary functions. As an example, we choose 



/ = aisn(a 1 x + C\t + d\, mi) + a 2 sn(a 2 x + c 2 t + d 2 , m 2 ), 

h = (3sn(by + d 3 ;m 3 ), (27) 

where b, (3, a*, a i; c, and dj are arbitrary constants (i = 1,2; j = 1,2,3) and 
rn^'s are modulus parameters. Then 

« = ^. (28) 

where, q± = 2[(aiaicn(«i; rai)dn(iii; mi) + a 2 a 2 cn(u 2 ; m 2 )dn(u 2 ; m 2 )) 
(3bcn(u 3 ; m 3 )dn(u 3 ; m 3 )], q 2 = (c+aisn^; mi)+a 2 sn(-u 2 ; m 2 )+f3sn(u 3 ; m 3 )) 2 , 
U\ = a\X + C\t + d\, u 2 = a 2 x + c 2 t + d 2 and u 3 = by + d 3 with corresponding 
expressions for v(x,y,t). We choose the constant c such that |c| > \a± + 
a 2 + f3\ for nonsingular solutions. The profile of the above solution for the 
parametric choice a\ = a 2 = (3=1, a\ =0.5, a 2 =0.5, b =0.4, c=-4, ci=-l, 
c 2 =-2, d\=d 2 =d 3 =0, m\= 0.2, m 2 = 0.3, m 3 = 0.4 at time t=0 is shown in 
Fig. 2a. As m 1? m 2 , m 3 — > 1, the above solution given by Eq. (1251) degenerates 
into a (2,1) dromion solution given by 



2(ctiaisech 2 Mi + a 2 a 2 sech 2 'u 2 )6sech 2 -U3 
(c + aitanhu! + a 2 tanhu 2 + /3tanhu 3 ) 2 



(|c| > |ai + a 2 + /9|) (29) 



where u\ = a\X + C\t + d\, u 2 = a 2 x + c 2 t + <i 2 and u 3 = by + d 3 . The solution 
is plotted for the parametric choice «i = a 2 = (3=1, a\ =0.5, a 2 = 0.8, 6 =0.4, 
c=-4, Ci=-1, c 2 =-2.5, c?i=(i 2 =(i3=0, for various values of t in Figs. 2(b)-2(d) 
in order to bring out the dromion interaction clearly. 



3.1.3 Asymptotic analysis for (2,1) dromion solution 

Since both the dromions in Figs. 2(b)-2(d) corresponding to the solution (1291) 
are travelling along the x-direction, following one another, it is enough to do 
this analysis for y = (and d 3 = so that u 3 = 0). Similar analysis holds good 
for any other value of y and d 3 ^ 0. This restriction corresponds to the cross 
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section of dromions, which are essentially solitons. We analyse the limits t — > 
— oo and t — > +00 separately so as to understand the interaction of dromions 
centered around u\ ~ or m 2 « 0. Without loss of generality, let us assume 
c\ > c 2 and ai < a 2 . Then, we find in the limit t — > ±00, Mi = aix + Cit + di 
and m 2 = a 2 x + c 2 t + d 2 take the following limiting values. 

1) As t -> -00: 

Mi ~ 0, M 2 — >■ +OO 
M 2 ~ 0, U\ — > — CXD 

2) As t -»• +00: 

Ml ~ 0, M 2 — > — OO 
M 2 ~ 0, Mi — > +OO 



1. Before interaction (as t — > —00): 

For Mi ~ 0, m 2 — » +00 and «i = a 2 = /3 = 1, the (2,1) dromion solution (J29 
becomes (soliton solution corresponding to dromion 1) 



2ai6 



-sech 2 (Mi + Si), Si 



'c + 2 



(30) 



For m 2 ~ 0,Mi 
dromion 2) 



c(c + 2) 

—00, ff29l) becomes (soliton solution corresponding to 
2a 2 6 



c(c-2) 



sech (m 2 + S 2 ), 5 2 



c-2 



(31) 



2. After interaction (as t — > +00): 

For Mi w 0,m 2 — > —00, fT29l) becomes (soliton solution corresponding to 
dromion 1) 



2aib 



-sech 2 (Mi + S2). 



(32) 



For m 2 w 0, Mi 
dromion 2) 



+00, 



c(c-2) 

becomes (soliton solution corresponding to 



2a 2 6 
c(c + 2) 



sech 2 (M 2 + Si). 



(33) 



The above results can be interpreted in the following way. Before interaction, 
the dromion 1 has a larger amplitude, travelling slower and the dromion 2 has 
a shorter amplitude, travelling faster. During interaction (at time t = 0), an 
exchange of energy between the dromions take place. This results in the gain 
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in amplitude of dromion 2 and fall in amplitude of dromion 1. But there is no 
change in velocity of the dromions and there is only a change in phase. 



3.1.4 Bounded multiple solitary waves 

One can also construct bounded two solitary waves by choosing 



/ = asn(ax + c\t + d\\ mi), 

h = Ptsn^y + d 2 ;m 2 ) + (3 2 sn(b 2 y + d 3 ; m 3 ), (34) 

where a, a, ci, /3j, bi, and dj (i = 1, 2; j = 1, 2, 3) are arbitrary constants and 
mj's are modulus parameters. Then 

? = -, (35) 

where q 3 = 2aacn(w 1 ; m 1 )dn(M 1 ; mi)[/3i6icn(w 2 ; m 2 )dn(-u 2 ; m 2 ) + j3 2 b 2 cia(u 3 ; m 3 ) 
dn(M 3 ;m 3 )], q A = (c + asn(«l; ml) + /3isn(n2; m2) + /5 2 sn(«3; m3)) 2 , Ui = 
aix + Cit + di, u 2 = biy + d 2 and u 3 = b 2 y + d 3 with corresponding expressions 
for v(x,y,t). Here \c\ > \a + fii + [3 2 \. The profile of the above solution for 
the parametric choice a = (3\ = (3 2 = 1, a =0.5, b\ =2.5, b 2 =1.7, c=-4, Ci=l, 
c 2 =-2, ^1=^2=^3=0, mi= 0.2, m 2 = 0.3, m 3 = 0.4, t=0 is shown in Fig.3(a). 
As mi, m 2) m 3 -> 1, the above solution, Eq. (1351) . degenerates into a bounded 
two dromion solution given by 

2a asech 2 ui(/9i6isech 2 u 2 + /3 2 fe 2 sech 2 w 3 ) 
(c + atanhiii + /3itanhu 2 + /3itanh%) 2 ' 

where u\ = ax + Cit + di, m 2 = 612/ + <^2 an d u 3 = b 2 y + d 3 . The above solution 
for the parametric choice a = fli = /3 2 = 1, a =0.5, 61 =2.5, 6 2 =1.7, c=- 
4, Ci=l, c 2 =-2, di=d 2 =d 3 =0, is shown in Fig.3(b). Here both the dromions 
travel with equal velocity along the x direction. Since they move parallel to 
each other there is no interaction between them. 



3.2 (N,M) dromion solution 



Proceeding in a similar way as above, one can generate (N,M) dromion solution 
by choosing 
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N 

/ = X) a i sn ( a i x + erf + da] rrii), 
1=1 

M 

^ = X (3isn(bjy + d 2j ; mj), 
3=1 



(37a) 
(37b) 



where a», a», 6j, Cj, d^, c?2j are arbitrary constants and m,, are modulus 
parameters and substituting them in the expression for q and v in (1221) . As 
the analysis is similar to the above, we do not present further details here. 



3.3 Singular solutions 



While presenting the explicit dromion solutions above, we have always as- 
sumed the magnitude of the parameters c in the solution (122]) to be greater 
than a constant value in order to have nonsingular solutions. For example, for 
the (1,1) dromion solution fl25|) . we have taken |c| > |ot + /3|. For |c| < |a + /3|, 
the solution (1251) becomes singular and it is illustrated in Figs.4(a)-4(c) for 
different instants of time. Here we have chosen a = /3 = 1, c=l, a — 0.5, b = 
0.4, c = 1, d± — di = 0, a\ = a.-} — 1- 



4 Solutions corresponding to product of arbitrary functions 



In this section, we will assume c 2 = c 3 = and c 4 = 1 in (1T51) or (f2"0"|) - (12T|) and 
call c\ = c. Then the solution can be written as 



-2cf x g y 

(fg + c 

2 fig 2 2f xx g , f t + f : 



2 , (38a) 



(f9 + c) 2 (fg + c) 3f 
Now we will look for periodic and localized solutions. 



+ - — f^. (38b) 



4-1 Harnessing of arbitrary functions and novel solutions of (2+1) dimen- 
sional K-dV equation 



As before, we choose the arbitrary functions / and g to be Jacobian elliptic 
functions, namely sn or cn functions, 

/ = asn(ax + c\t + d\\ mi), g = /3sn(by + d 2 ; m 2 ) (39) 
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so that 

q(x,y,t) = ^, (40) 
Qe 

where 05 = — 2 cot (3 a bcn.{u\, mi)dn(«x; mi) cn(it 2 ; m 2 )dn(it 2 ; m 2 ), qe = (c + 
a/3sn (mi; mi)sn(u 2 ; m 2 )) 2 , «i = ax + C\t + di and u 2 = by + d 2 . In Eqs. f l39l 
and PU|) . the quantities mx and m 2 are the moduli of the respective Jaco- 
bian elliptic functions while a, (3, a, 6, c, ci, di and d 2 are arbitrary constants 
with |c| > The profile of the above solution for the parametric choice 

a — (3 — 1, a=0.5, 6=0.4, c=-4, Ci=-1, di=d 2 =0, 0.2, m 2 = 0.3 at time 
t=0 is shown in Fig. 5(a). It can be observed that here again the periodic wave 
moves along the x-direction only. 



4-1.1 (1,1) dromion solution 

As a limiting case of the periodic solution given by Eq. (l4"0l) when mi, m 2 
— > 1, the above solution degenerates into an exponentially localized solution 
(dromion). The limiting form corresponding to the (1,1) dromion take the 
expression 

—2c a (3 absech 2 (ax + cit + d 1 )sech 2 (6y + d 2 ) 



(c + a (3 tanh(ax + C\t + di)tanh(by + d 2 )) 2 
with |c| > | a/? | and the variable v then takes the form (using expression (I38bl )) 



2a 2 a 2 (3 2 sech 4 (ax + c x t + aMtanh 5 \by + d 2 ) 



(c + a /3 tanh(ax + c\t + di)tanh(6y + d 2 )) 2 
^Aa 2 a (3 sech 2 (ax + cit + di)tanh(ax + Cit + di)tanh(6y + d 2 ) 
(c + a /? tanh(ax + Cit + c?i)tanh(6?y + d 2 )) 

— [c x + 4a 3 tanh 2 (ax + c x t + di) - 2a 3 sech 2 (ax + c\t + di)]. (42) 

Schematic form of the (1,1) dromion for the parametric choice a = (3 = 1, 
a— 0.5, 6=0.4, c=-4, Ci=-1, d!=d 2 =0 at time t=0 is shown in Fig. 5(b). Again 
note that the dromion travels along the x-direction. 



4-1.2 More general periodic solution and (2,1) dromion 

Next, we obtain more general periodic solution by choosing further general 
forms for the arbitrary functions. As an example, we choose 



/ = a\sa{a\X + c\t + d\\ m\) + a 2 sn(a 2 x + c 2 t + d 2 ; m 2 ), 

g = f3sn(by + d 3 ;m 3 ), (43) 
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where b, (3, a>i, Ci and dj are arbitrary constants (i = 1,2; j = 1,2,3) and 
m^s are modulus parameters. Then 

, = * (44) 
<Zio 

where q 9 = — 2c (3 [(aia 1 cn(w 1 ; mi)dn(wi; mi) + a 2 a 2 cn(-u 2 ; m 2 )dn(-u 2 ; m 2 )) 
6cn(u 3 ; m 3 )dn(w 3 ; m 3 )], q w = [c + /3(aisn(ii 1 ; mi) + a 2 sn(w 2 ; m 2 ))sn(%; m 3 )] 2 , 
= + c\t + di, u 2 = a 2 x + c 2 £ + d 2 and u 3 = by + d 3 with corresponding 
expressions for v(x,y,t). Here |c| > |(ai + a 2 )(3\. The profile of the above 
solution for the parametric choice a\ = a 2 = (3 — 1, a\ =0.5, a 2 =0.8, b =0.4, 
c=-4, ci=-l, c 2 =-2.5, c?i=(i 2 =(i 3 =0, mi= 0.2, m 2 = 0.3, m 3 = 0.4 at time 
t=0 is shown in Fig. 6 (a). As mi, m 2 , m 3 — > 1, the above solution Eq. (1441) 
degenerates into a (2,1) dromion solution given by 

—2c 6 /3(aiaisech 2 Mi + a 2 a 2 sech 2 w 2 ) sech 2 w 3 



(c + /3(aitanhui + a 2 tanhu 2 )tanhu 



3 ) 



\c\ > |(ai+a 2 )/?|, where Mi = aiX+Cit+di, u 2 = a 2 x+c 2 t+d 2 and u 3 = by+d 3 . 
The dromion interaction for the parametric choice a± — a 2 — (3 — 1, a,\ =a 2 
=0.5,6 =0.4, c=-4, Ci=-1, c 2 =-2, rfi=(i 2 =c? 3 =0, is shown in Figs.6(b)-6(d) for 
different time intervals. 



4-1.3 Asymptotic analysis for (2,1) dromion solution corresponding to Eq. jj^5\) 

Proceeding with the analysis as was done in the previous section, we find that 
both the dromions are of different amplitudes and travelling with different 
velocities along the x-direction. At time t = 0, there is an interaction between 
the dromions. After interaction, we find that there is no change in amplitude 
or velocity of the dromions . 

4-1-4 Bounded multiple solitary waves 

One can also construct bounded two solitary waves as before by choosing 



/ = asn(ax + c±t + d\, mi), 

g = PiSiifay + d 2 ;m 2 ) + f3 2 sn(b 2 y + d 3 ; m 3 ), (46) 

where a, a, Ci, 6j and dj are arbitrary constants (i — 1, 2; j — 1, 2, 3) and m^'s 
are modulus parameters. Then 

q = ^, (47) 

<?12 

where qu = —2c a a cn(«i; mi)dn(«i; mi)[/?i6icn('U 2 ; m 2 )dn(w 2 ; m 2 )+/5 2 6 2 cn('U 3 ; 
m 3 )dn(w 3 ;m 3 )], q u = (c + asn(«i; mi)[(3 x sn(u 2 \ m 2 ) + (3 2 sn(u 3 ; m 3 )]) 2 , U\ = 



13 



a\x + c\t + d\, u 2 = b±y + d 2 and us = b 2 y + d 3 with corresponding expressions 
for v(x,y,t). Again here |c| > \a(fti + (3 2 )\ for nonsingular solutions. The pro- 
file of the above solution for the parametric choice a = (3\ = 0% = 1, a =0.5, 
bi =2.5, b 2 =1.7, c=-4, Ci=l, c 2 =-2, di=d 2 =d 3 =0, 0.2, m 2 = 0.3, m 3 = 
0.4, t=0 is shown in Fig. 8a. As mi, m 2 , m 3 — > 1, the above solution given by 
Eq. (|47p degenerates into a bounded two dromion solution given by 

— 2caasech 2 -u 1 (/? 1 & 1 sech 2 w 2 + /3 2 fr 2 sech 2 -u 3 ) 
(c + atanhui(,5itanhu 2 + /3 2 tanlru 3 )) 2 

|c| > + f3 2 )\, where u\ = ax + c±t + di, u 2 = b±y + d 2 and U3 = b 2 y + d% 
with the parametric choice being a = fli = /3 2 = 1, a =0.5, b\ =2.5, b 2 =1.7, 
c=-4, Ci=l, c 2 =-2, rf!=(i 2 =(i3=0. The (2,1) dromion is shown in Fig. 7(b). The 
two dromions evolve quite similar to the pattern of Fig. 3(b) in the x-direction. 



5 Generalized Nizhnik-Novikov-Veselov (NNV) Equation and con- 
struction of solutions 

The generalized Nizhnik-Novikov-Veselov (NNV) equation is a symmetric gen- 
eralization of the K-dV equation in (2+1) dimensions and is given by 



0, (49a) 
v v , (49b) 
w x . (49c) 

Here a, b, c and d are parameters. This equation which is also known to 
be completely integrable has been investigated and exponentially localized 
solutions have been generated [4"]I8"] . We now apply the Painleve truncation 
approach to obtain more general solutions. For this purpose, we again truncate 
the Laurent series at the constant level term to get the transformation 



u = uq4> 2 + U\(j) 1 + u 2 , (50a) 

v = v (p- 2 + v^' 1 + v 2 , (50b) 

w = wq4>~ 2 + Wicj)^ 1 + w 2 . (50c) 

Again considering the vacuum solutions of the form 

u 2 = 0,v 2 = v 2 (x,t),w 2 = w 2 (y,t), (51) 



where v 2 (x,t) and w 2 (y,t) are arbitrary functions in the indicated variables, 
and collecting the coefficients of different powers of </> as before and solving 
the resultant equations, we obtain the following results: 



u t + au xxx + buyyy + cu x + duy — Sat ] x u — 3avu x 

—3bw y u — 3bwu y = 
u x = 

Uy = 
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M = 2<f> x <j)y, ui = -2(p xy , (52a) 

v = 2<j>l, v x = -2(j) xx , (52b) 

w = 2(f>y, wi = -2<p yy . (52c) 

Also v 2 and w 2 are related through </> by the relation 



^ _ <j>t + + {$>x$xxy — fyxxfyxy) b{(j) y (j) X yy — (j>yy(f>xy) 

3a<j) x 4> x 4>y acj) 2 x 

+ c_ + dcp^ _ bcf)yW2 ^ 
3a 3acj) x acp x 

Making use of the expression (1521) and (I53I) in (1501) and using the resultant 
forms in the generalized NNV equation (1491) . we obtain the following system 
of equations which is triliner in cj) (with the coefficient w 2 also present), 



( t ) xx{ ( t ) y ( t ) xyy ~ (ftyytpxy) + 'Pxi'Pyy^xxy ~ ( f ) y < fixxyy) — 0, (54a) 
My{4>yt + CKpyxxx + Hyyyy) ~ 4>y{4>t + a<j) xxx + b<t> yyy ) (j) xy 

3d{^) X ^) XX y <fixx4 > Xy)4 > Xy 4" 3b<fiy((J)y(j) X yyy ^yyy^Xy) 

yy4*xyy ^y^xxyy 

) + d(f>y(4> x (f)yy ~ <f>y(fi X y) 

-3b(j) y [(j) x ((j)yyW2 + 4>yW 2 y) - (f)y(j) X yW2} = 0, (54b) 

w 2x = 0. (54c) 

One can observe that the above set of equations ( 1541) admits a more gen- 
eral form of <p involving two arbitrary functions f(x,t) and g(y,t), and four 
arbitrary parameters, ci, c 2 , c 3 and c 4 as 

(f)(x, y, t) = a + c 2 f(x, t) + c 3 g(y, t) + c 4 f(x, t)g(y, t). (55) 

Here v 2 and w 2 take the forms (from (l53l) ) 



ft + af xxx + cf x + (eg + c 4 f) , s 

V2 = 5—? , (56a) 

9t + bg yyy + dg y - (pa + c 4 #) 

^2 = ^7 • (56b) 

3bg y 

Here ci, C2, C3, and C4 are arbitrary parameters. Note here that unlike the case 
of (2+1) dimensional KdV equation only two arbitrary functions are allowed 

m am. 



Now substituting the above form (1551) for the manifold <j)(x, y, t) into the trun- 
cated Painleve series (!50|) for the functions u(x,y,t), v(x,y,t) and w(x,y,t), 
along with the expressions for the coefficient functions u , v , w Q , u±, vi, 101, 
v 2 and w 2 given above, we finally obtain the solution to Eqs. (l4"9l as 
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.,( T t\ = 2(c 2 c 3 - c 4 c x )f x g y 

1 ' V ' } h + c 2 f(x, t) + c 3 g(y, t) + c,f(x, t)g(y, t)] a ' 1 J 



v(x,y,t) 



w(x,y,t) 



(c 2 c 3 - c 4 Ci) 7^ 0, 

2(c 2 + c 4 ^) 2 /, 2 



[c\ + c 2 f(x, t) + c 3 g(y, t) + c 4 f(x, t)g(y, t)} 2 

2(c 2 + c 4 g)f xx 

ci + c 2 /(x, £) + c 3 g(y, t) + c 4 /(x, t)g(y, t) 

. ft + fxxx + Cf x + (C 3 + C 4 /) 



2(c 3 + c 4 /)V 



[ci + c 2 f(x, t) + c 3 g(y, t) + c 4 f(x, t)g(y, t)} 2 

2(c 3 + c 4 f)g yy 

ci + c 2 f(x, t) + c 3 g(y, t) + c 4 /(x, t)g(y, t) 

36^ 



(57b) 



(57c) 



One can check that the above solutions of NNV equation reduces to the 
solutions of (2+1) dimensional KdV equation for the parametric choice b = 
c = d = and also by considering the truncated w series to be zero. One 
may note that the above form of solution coincides with the universal form of 
solutions reported by Tang et al. [5] and the special case c\ — A (A is arbitrary 
constant), c 2 = c 3 = was studied by Peng [8]. One can obtain periodic and 
localized dromion solutions here also following the procedure discussed in the 
earlier sections for the (2+1) dimensional KdV equation. 



6 Discussion 



In this paper, we have investigated the (2+1) dimensional K-dV equation and 
obtained a four parameter solution involving three arbitrary functions by us- 
ing the Painleve truncation approach. For different choice of parameters, we 
have generated two broad classes of localized coherent structures and elliptic 
function periodic wave solutions. In particular, we have shown that the equa- 
tion is trilinearizable. Even more general solutions can be generated from the 
trilinear equations which remains to be investigated. The existence of lower 
dimensional arbitrary functions helps us to construct novel localized solutions 
and study their interactions. We have also extended the approach to general- 
ized NNV equation and pointed out how similar solutions as that of the (2+1) 
dimensional KdV equation can be obtained. 
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A One dromion solution through Hirota Bilinearization 



Here we briefly point out how the (1,1) dromion solution can be obtained 
through Hirota bilinearization method [I]. To bilinearize equation ([3]), we make 
the transformation 

u = -2d xy (\og(j)), v = -2d xx (\og(j)), (A.l) 

which can be identified from the Painleve analysis. The resultant bilinear form 
is given by 

(D y D t + DlD y )(j) -0 = 0, (A26) 

where D's are the usual Hirota operators. To generate a (1,1) dromion, one 
considers the ansatz 

(j) = l + e^ 1 + e^ 2 + Ke^ 1+1p2 , (A3) 

where 

-0 1 = kix - k\t + 5 U (A4a) 

^2 = hy + S 2 - (AAb) 

Here ki, h, 5%, 82 and K are arbitrary constants. The (1,1) dromion is given 

by m 

2fc 1 / 1 (l-AT)e^+^ 
U (1 + + + Ke*i+4»)* ' 1 } 

This is a special case of the dromion solution we have obtained in (123]) with the 
constants a = (3 = I, a = y, & = |, c = -75=5*, ci = -y, <2i = |log^ + y, 
a? 2 = |log— + y. Thus the standard dromion solutions become special cases 
of the general solutions obtained in section 3. 

Figure captions 

Figure 1 (a) Elliptic function solution (124"]) for q(x, y, t),(b) Localized dromion 
solution (1231) for q(x, y, t) and (c) the solution corresponding to (1231) for v(x, y, t) 
Figure 2 (a) Elliptic function solution ([23]) for q(x,y,t), (b-d) (2,1) dromion 
solution (1231) for g(x, y, t) and its interaction at time units (b) t = —4, (c) 
t = and (d) t = 4 

Figure 3 (a) Elliptic function solution (1331) . (b) bounded two dromion solu- 
tion flUD 

Figure 4 Singular solution for (1231) at different instants of time (a)£=-l, 
(b)t=0, (c)t=l 

Figure 5 (a) Elliptic function solution (T431) for g(x, y, t) (b) Localized dromion 
solution (|4T|) for q(x,y,t) (c) the solution corresponding to (142]) for v(x,y,t) 
Figure 6 (a) Elliptic function solution (144]) for q(x,y,t), (b-d) (2,1) dromion 
solution (1431) for q(x,y,t) and its interaction at time intervals (b) t = —4, (c) 
t = and (d) t = 4 
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Figure 7 (a) Elliptic function solution (H7I) . (b) bounded two dromion solu- 
tion flUD 
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